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■ We present a method to calculate vertical profiles of particle size distributions in 

. condensation clouds of giant planets and brown dwarfs. The method assumes a bal- 

ance between turbulent diffusion and sedimentation in horizontally uniform cloud decks. 
| Calculations for the Jovian ammonia cloud are compared with results from previous 

methods. An adjustable parameter describing the efficiency of sedimentation allows the 
; new model to span the range of predictions made by previous models. Calculations 

for the Jovian ammonia cloud are consistent with observations. Example calculations 
are provided for water, silicate, and iron clouds on brown dwarfs and on a cool extra- 
solar giant planet. We find that precipitating cloud decks naturally account for the 
characteristic trends seen in the spectra of L- and T-type ultracool dwarfs. 
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Introduction 



The visual appearance and spectrum of every solar system body with an atmosphere depends 
strongly upon the character and distribution of atmospheric condensates. This is particularly true 
for the giant planets where optically thick cloud decks dominate the appearance of the planets at 
most continuum wavelengths in both the reflected solar and the thermal infrared. Condensates also 
play a role in controlling the spectra of at least some brown dwarfs and most extrasolar giant planets. 
Indeed one suggested classification scheme (Sudarsky, Burrows & Pinto 2000) for extrasolar planets 
hinges on the specific atmospheric condensates present. Yet despite the importance of condensates 
there exists no simple model for predicting the parameters most relevant to radiative transfer: the 
vertical profile of condensate mass and its distribution over particle size. 
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Chemical equilibrium models (e.g., Lewis 1969; Fegley & Lodders 1994) predict which species 
are expected to condense in an atmosphere, yet they provide no guidance as to the expected particle 
sizes. Other models (e.g., Rossow 1978; Lunine et al. 1989; Carlson, Rossow, Sz Orton 1988) predict 
some parameters, but lack a simple, self-consistent recipe for exploring the possible phase space in 
which clouds might exist. 

A single example motivates the need for cloud models in substellar atmospheres. With increas- 
ingly later spectral type the warm L-dwarfs become progressively redder in their J — K color (e.g., 
Kirkpatrick et al. 1999; Martin et al. 1999; Fan et al. 2000). Spectral fitting and models (e.g., 
Leggett, Allard, k Hauschildt 1998; Chabrier et al. 2000; Marley 2000) demonstrate that this is 
due to the progressive appearance of more silicate dust in the cooling brown dwarf atmospheres. 
Yet the cooler T-type brown dwarfs, like Gliese 229B, have blue colors in J — K (e.g., Leggett et 
al. 1999; Tsvetanov et al. 2000). The spectra and colors of these cool brown dwarfs can only be 
fit by atmosphere models that assume the silicate dust has settled below the visible atmosphere 
(e.g., Allard et al. 1996; Marley et al. 1996; Tsuji et al. 1996). Models in which the dust does not 
settle (Chabrier et al. 2000) produce T-dwarf colors that are at least 2 to 3 magnitudes redder than 
observed. Marley (2000) has demonstrated that a simple model in which all clouds are a single 
scale-height thick can explain this behavior, but the assumed distribution was prescribed rather 
than being calculated from any model physics. Correct modeling of the atmospheres of cooling 
brown dwarfs and the ultimate assignment of an effective temperature scale to the L-dwarf spectral 
sequence (e.g., Kirkpatrick et al. 1999; Basri et al. 2000) requires a characterization of clouds. The 
ideal model would have a small number of free parameters, predict the vertical distribution and 
particle sizes of the condensates, and yet be simple enough to be included into model atmosphere 
codes that iteratively search for self-consistent atmospheric structures. No such ideal model yet 
exists. 

We aim to fill this void by presenting a simple model describing precipitating clouds in sub- 
stellar atmospheres. We limit our treatment to condensation clouds, and hence do not consider 
photochemically driven hazes likely to appear in illuminated stratospheres. We depart from previ- 
ous work by explicitly treating the downward transport of raindrops with sizes greater than that 
predicted from the convective velocity scale. Including rainfall produces clouds of thinner vertical 
extent, which can better reproduce observations of Jupiter's ammonia cloud. The resulting model 
is general enough to be applied to iron and silicate clouds appearing in brown dwarf atmospheres 
(e.g., objects with effective temperatures, T e g- ~ 1500 K) as well as the atmospheres of cool extra- 
solar giant planets (T e g- ~ 400 K) in which water clouds dominate the atmosphere. The few free 
parameters in the model can produce clouds with dramatically different characteristics; ultimately 
observations will constrain these parameters and hopefully provide information on the underlying 
atmospheric dynamics and cloud physics. 

In this paper we first summarize previous cloud modeling efforts, then describe the new model. 
We use the ammonia cloud of Jupiter as a framework for describing the model physics and evaluating 
the model performance. Finally we illustrate model applications by considering water, silicate, and 
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iron clouds in the atmospheres of brown dwarfs and a cool extrasolar giant planet. 

2. Previous Models 

A great range of models have been used to represent clouds in the terrestrial atmosphere, which 
vary in the complexity by which atmospheric dynamics and cloud microphysics are treated. The 
most detailed models simulate three-dimensional cloud-scale motions and resolve the size distribu- 
tions of cloud droplets (and the aerosols on which they form) and treat the interactions between 
dynamics, microphysics, and radiative transfer. The computational demands of such complex 
models limit their domain sizes to a few kilometers in each dimension. Present global-scale (gen- 
eral circulation) models greatly simplify the representation of clouds by parameterizing cloud-scale 
motions as well as cloud microphysical processes, and such simplifications lead to profound uncer- 
tainties in climate predictions from their simulations. Both types of models, as well as a range 
of intermediate models, can be considered appropriate for modeling the terrestrial atmosphere by 
virtue of the wealth of observational data available to constrain them; whether or not the unknowns 
in such models are uniquely constrained by the data constitutes a debate beyond the scope of this 
study. 

The relative scarcity of observational data for clouds in extraterrestrial atmospheres is far less 
constraining. Leading uncertainties include the characteristics of atmospheric dynamics and the 
populations of nuclei upon which cloud droplets form. Hence, we consider it appropriate to model 
extraterrestrial clouds using much simpler treatments. 

Perhaps the simplest approach to modeling clouds is through a Lagrangian parcel model, in 
which the base of a cloud appears where the adiabatic cooling of an air parcel in an updraft results 
in saturation (ignoring any supersaturation associated with barriers to cloud droplet formation). 
Further cooling condenses vapor in excess of saturation onto cloud particles. The particles grow 
through condensation and coalescence until their sedimentation velocities exceed the updraft speed, 
and then fall out of the parcel. A number of problems arise in the formulation of updraft parcel 
models, among them: ignoring parcels in downdrafts, treating the mixing between parcels, treating 
the source of condensates into a parcel due to sedimentation from above, and determining updraft 
speeds. 

Another simple approach, which we employ here, is through a one-dimensional Eulerian frame- 
work, in which turbulent diffusion mixes a condensable vapor upwards, while maintaining a constant 
mixing ratio (equivalently, mole fraction) below the cloud. Temperature and hence the saturation 
mixing ratio in the air column decrease with altitude, and the cloud base again appears where 
the saturation mixing ratio matches the sub-cloud mixing ratio. Above the cloud base, turbulent 
diffusion works toward maintaining a constant total mixing ratio (q t = q v + Qc), which is the sum of 
the vapor (q v = moles of vapor per mole of atmosphere) and condensate (q c = moles of condensate 
per mole of atmosphere) mixing ratios, while sedimentation reduces qt by transporting condensate 
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downward. Note that by ignoring horizontal variability, any differences between (cloudy) updrafts 
and (potentially cloud-free) downdrafts are neglected. 

A number of models for tropospheric condensation clouds have appeared in the planetary and 
astrophysical literature. Here we summarize a selection of them that contribute to the present 
work. 



Lewis (1969) represents a foundation in the study of tropospheric clouds in the giant planets. 
In that work the term "precipitation" is used in the narrow sense used by chemists, in which 
condensates appear where the local saturation vapor pressure is exceeded by the actual vapor 
pressure, rather than in the broader sense used by meteorologists, which additionally denotes 
sedimentation of the condensates (hereafter we use the term in this broader sense). Although 
there is no mention of sedimentation by Lewis, the treatment does imply certain assumptions. 
Starting below the cloud base and working upwards, at each computational level the Lewis model 
assumes that all the condensate remains at the level where it appears. Considered in the framework 
of a parcel in an updraft, the Lewis model assumes that all condensate rains out with a fallspeed 
matching the updraft velocity. Were sedimentation slower, condensate would be transported upward 
(as discussed by Weidenschilling Sz Lewis 1973); were sedimentation faster, condensate would be 
transported downward. Hence, the Lewis (1969) assumption regarding sedimentation is an unstated 
compromise between those two possibilities. 

We implement the Lewis model by starting below the cloud base (where q c = and q v = % e iow) 
and condensing all vapor in excess of saturation at each successive level upward: 



where z is altitude and q s is the vapor mole fraction corresponding to the saturation vapor pressure 
at that altitude. The first and second cases on the right hand side correspond to cloud-free and 
cloudy conditions, respectively. Note that under all conditions qt(z) = q v (z — Az) in the Lewis 
model, reflecting the assumption that only vapor is transported upwards. 

Beyond this simple model, Lewis (1969) considered the partitioning of chemical species in 
some detail, and also calculated pseudo-adiabatic lapse rates. Here we simply assume that each 
condensate results from the saturation of a single condensable, and fix the lapse rate as input from 
observations or an external model. 

For an example, we calculate an ammonia cloud profile from the Lewis model (Figure 1) using 
the Jovian temperature profile from Voyager (Lindal et al. 1981), the relation for vapor pressure 
given in Appendix A, and a sub-cloud mole fraction of 3 x 10~ 5 (a wide range of abundances below 
the expected base of the Jovian ammonia cloud have been reported; we adopt the value at 0.6 bars 



2.1. Lewis (1969) 



q c (z) = max(0, q v (z - Az) - q s {z)) 
q v {z) = min(q v (z - Az),q s (z)) 



(1) 
(2) 
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retrieved by Kunde et al. 1982 for the Northern Equatorial Belt, which also agrees with the best-fit 
values of Carlson, Lacis, & Rossow 1993 and Brooke et al. 1998). The cloud base appears at 0.42 
bars, where the temperature is 129 K. Although absent in the figures of Lewis (1969) (likely due to 
reduced vertical resolution), in our interpretation of that model the vapor is not entirely depleted 
in the lowest reaches of the cloud (where q c < q t ), hence q c increases above the cloud base. Such 
an increase is found in terrestrial clouds of moderate vertical extent, where q c < q t , and hence q c 
increases with altitude throughout their depths. However, at greater altitudes in this deep ammonia 
cloud, the vapor is so effectively depleted by condensation at the low temperatures that q v <C qt, 
leading to a cold degeneracy: q c (z) ~ qt(z) = q s (z — Az), in which decreasing temperatures result 
in q c diminishing with altitude. Note that the condensate abundance drops off rapidly above ~0.13 
bars due to increasing temperatures. Hence, the temperature minimum quite reasonably produces 
a cold-trap in the Lewis model. 

Condensate particle sizes, the other ingredient needed for predicting cloud opacity, are not 
considered by Lewis (1969) or Weidenschilling & Lewis (1973). 

2.2. Carlson et al. (1988) 

In their theoretical characterization of cloud microphysics of the giant planets, Carlson et al. 
(1988) employ the formalism of Rossow (1978) to calculate time constants for droplet condensation 
within cloudy updrafts (assuming a supersaturation of 10~ 3 ), droplet coalescence (assuming the 
mean collision rate is described by particles with a mass ratio of 2), and sedimentation through 
an atmospheric scale height. From these time constants, estimates are made of the predominant 
size of cloud particles at cloud base for a number of condensates. For the Jovian ammonia cloud, 
Carlson et al. estimate a mass-weighted droplet radius of ~ 10 — 30 fim. 

Carlson et al. (1988) make no attempt to calculate vertical profiles of condensate mass. For 
profiles of vapors that condense into multiple forms (such as ammonia, which can also condense 
onto a cloud of NH4SH below the ammonia cloud), saturation is assumed above the cloud base. 

A shortcoming to the approach of Carlson et al. (1988) is that their microphysical time 
constants strongly depend on a number of uncertain factors, chief among them completely unknown 
supersaturations, which govern droplet growth rates due to condensation. Supersaturations in a 
cloudy updraft are determined by balance between the source due to adiabatic cooling, and the 
sink due to condensation. Uncertainties in updraft speeds and the populations of condensation 
nuclei (and hence cloud droplets) both contribute to the uncertainty in supersaturations realized in 
extraterrestrial clouds. Furthermore, the time constants Carlson et al. (1988) use for gravitational 
coalescence assume that the collection efficiency is unity, and those for sedimentation effectively 
assume a fixed width of the size distribution. Rather than attempting to constrain the many degrees 
of freedom using such a detailed approach, we choose instead to reduce the number of assumptions 
by simplifying the description of cloud microphysics. 
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2.3. Lunine et al. (1989) 

Lunine et al. (1989) consider a range of possible iron and silicate clouds in brown dwarfs; the 
possibilities differ in the nature of the balance between sedimentation and turbulent mixing. The 
framework is based on a theoretical investigation into iron clouds deep in the Jovian atmosphere by 
Prinn & Olaguer (1981), which in turn draws on an analysis of sulfuric acid clouds on Venus (Prinn 
1974). These models represent a fleshing out of the discussion of vertical transport of condensates 
by Weidenschilling & Lewis (1973). 

Two fundamental cloud types are treated by Lunine et al. (1989). The first is "dust-like" (using 
the terminology of Prinn & Olaguer 1981) in which cloud particles grow and efficiently sediment 
out, resulting in relatively thin clouds limited by the local vapor pressure, as in the model of Lewis 
(1969). These dust-like clouds are assumed to prevail in the radiative region (stratosphere), where 
the temperature profile is stable and convection is suppressed. 

The second fundamental type in the Lunine et al. (1989) study is a tropospheric cloud, in 
which downward transport by sedimentation is opposed by upward transport due to turbulent 
mixing. For this cloud type, two variations are considered by Lunine et al. (1989). For the first 
variation, described as "frozen-in," cloud particles are so small that sedimentation is overwhelmed 
by upward transport due to turbulent mixing. In this case, the atmosphere is well-mixed with 
respect to condensate, and hence q c is independent of altitude above the cloud base. For the second 
variation, which is intermediate to the dust-like and frozen-in cases, particles grow large enough 
in "convective" clouds to develop appreciable sedimentation velocities, and their downward sedi- 
mentation is balanced by their turbulent transport upward. For their calculations of specific brown 
dwarf models, Lunine et al. (1989) consider only the two endmembers of their cloud spectrum, 
corresponding to dust-like and frozen-in clouds. 

Their intermediate case serves as a starting point for our model of condensate mass profiles. 
Our interpretation of the convective cloud model of Lunine et al. (1989) as applied to the Jovian 
ammonia cloud is shown in Figure 1. Note that we have refined that model slightly, allowing the 
atmospheric properties to vary with height above the cloud base, and relaxing their assumption 
that q c = q t . The condensate mass is seen to be significantly enhanced above the cloud base for that 
model: at the tropopause (where there is no cold-trap in this case) q c is enhanced a thousand-fold 
over that computed by the Lewis (1969) model. Thus the treatment in which Lunine et al. (1989) 
assume particle sedimentation to balance turbulent transport results in a cloud not so different 
from their frozen-in case (as depicted by the curve in Figure 1 labeled / ra i n = 0). Evidently the 
sedimentation in this convective cloud model is far less effective than that assumed by Lewis (1969). 
As described below, for our calculations of condensate mass profiles we extend the Lunine et al. 
(1989) approach by applying a scale factor to the particle sedimentation. 

For radiative calculations, Lunine et al. (1989) assume all particles in the frozen-in and dust- 
like clouds are 1 and 10 /j,m in radius, respectively. 
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2.4. Marley et al. (1999) 



The Marley et al. (1999) model of water and silicate clouds in extrasolar giant planets rep- 
resents a variation on the Lewis (1969) model. As in the Lewis model, the calculation of vapor 
pressure (or equivalently, q v ) assumes that any supersaturation is quenched locally by condensation 
(Equation 2). However, the calculation of the condensate mole fraction (q c ) represents a departure: 
instead of calculating it from the vapor pressure in the underlying layer from Equation 1, Marley 
et al. scale it to the local saturation vapor pressure with the following assumption: 



The parameter / s corresponds to the potential supersaturation prior to condensation. Marley et 
al. (1999) treat / s as an adjustable parameter, ranging from a baseline value of 0.01 to an extreme 
value of 1. The baseline model as applied to the Jovian ammonia cloud is shown in Figure 1. The 
condensate mass is seen to be diminished by a factor of ~100 relative to the Lewis (1969) model. 

Increasing / s to 1 results in a hundredfold enhancement of q c throughout the cloud compared 
to the baseline case, as shown in Figure 2a. The principle difference between that extreme and the 
Lewis (1969) condensate model is that for the former there is no regime near the cloud base akin 
to shallow terrestrial clouds, in which q c increases with altitude. This difference is attributable to 
a discontinuity of qt in the treatment of Marley et al.: below the cloud base qt{z) = q v (z — Az) as 
in the Lewis (1969) model, but above the cloud base qt(z) = (1 + fs) q s (z). 

For their calculations of cloud particle sizes, which are decoupled from their calculation of 
condensate mass, Marley et al. (1999) apply the formalism of Rossow (1978) to two atmospheric 
endmembers: first, a quiescent atmosphere, in which the mean particle size is determined from the 
condition that the sedimentation rate matches the faster of coagulation and condensation (at an 
assumed supersaturation of 0.01); and second, a turbulent atmosphere in which mixing is balanced 
by sedimentation. 

The first endmember is subject to a similar catalog of unconstrained assumptions as required by 
the treatment of Carlson et al. (1988), the most notable among them being the great uncertainty in 
the supersaturation driving droplet condensation. Also, the model physics underlying this quiescent 
atmosphere seems to be self-contradictory, on the one hand explicitly assuming that there is too 
little convection to regulate the maximum size of the droplets, yet on the other hand implicitly 
assuming that there is enough convection to supply the vapor necessary to drive condensational 
growth. 

However, the second case of Marley et al. (1999) requires significantly fewer assumptions, and 
is also appropriate to tropospheric condensation clouds. This second case serves as a starting point 
for the calculation of the cloud particle sizes in our model. 




(3) 
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3. The Present Model 

We model all condensation clouds as horizontally homogeneous (globally averaged) structures 
whose vertical extent is governed by a balance between the upward turbulent mixing of condensate 
and vapor (q t = q c + q v ) and the downward transport of condensate due to sedimentation: 

~ K lh ~ f rainW * Qc = ( 4 ) 

where K is the vertical eddy diffusion coefficient and / ra in is a new parameter that we have intro- 
duced, defined as the ratio of the mass- weighted droplet sedimentation velocity to w*, the convective 
velocity scale. We solve Equation 4 for each condensate independently, and hence ignore any mi- 
crophysical interactions between clouds. Equation 4 is an extension of Lunine et al.'s convective 
cloud model, relaxing their implicit assumptions / ra in = 1 and q c = Qt- 

The product f ra i n w* represents an average sedimentation velocity for the condensate, which 
offsets turbulent mixing and thereby leads to qt decreasing with altitude. The extreme case with 
no sedimentation to offset turbulent mixing (/ ra i n = 0) is equivalent the frozen-in endmember of 
Lunine et al. (1989) and the "dusty" models of the Lyon group (e.g., Chabrier et al. 2000). In this 
case the solution to Equation 4 is a well- mixed atmosphere (qt independent of altitude), which is 
seen in Figure 1 to loft even more condensate than the convective cloud of Lunine et al. (1989). 

We adopt / ra i n as an adjustable input parameter, which together with q c constrains the droplet 
size distributions. First we describe our calculation of q c , then the size distributions. 



3.1. Condensate Mass Profiles 

The eddy diffusion coefficient (K) for q t is assumed to be the same as that for heat as derived 
for free convection (Gierasch and Conrath, 1985): 



3 \H J \WaCp 

where the atmospheric scale height is given by H = RT/fig (for Jupiter we use g = 25 m s~ 2 ), L 
is the turbulent mixing length, R the universal gas constant, n the atmospheric molecular weight 
(2.2 g mol -1 assumed here), p a the atmospheric density, and c p the specific heat of the atmosphere 
at constant pressure (ideal gas assumed). Here we assume all the interior heat to be transported 
through the convective heat flux: F = <tT^ s , where a is the Stefan-Boltzmann constant and the 
effective temperature for Jupiter is T e g- = 124 K. In the general case, a profile of F is specified by 
an external model, which partitions the transport of interior heat between radiative and convective 
fluxes. The convective heat flux can be reduced further by other heat fluxes, such as para hydrogen 
conversion or latent heat release, as discussed by Gierasch and Conrath (1985). Beyond any uncer- 
tainty in the convective heat flux, the constant coefficient scaling the eddy diffusion coefficient is 
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only loosely constrained by observations. For our baseline model we use a coefficient of 1/3 based on 
previous modeling studies of the Jovian atmosphere (D. M. Hunten 1996, private communication), 
and consider the sensitivity of our model results to its value in a subsequent section. 

For freely convecting atmospheres, the mixing length is typically assumed to be the pressure 
scale height. However, in stable atmospheric regions, the mixing length will be diminished. We 
account for this reduction by scaling the mixing length to the local stability: 



where T and r a di a b are the local and dry adiabatic lapse rates, respectively, and A is the mini- 
mum scaling applied to L (we assume a value of 0.1). In the general case, convective heat fluxes 
are diminished in radiative regions; hence, we also assume an eddy diffusion coefficient no less 
than a prescribed minimum value (i^min = 10 5 cm 2 s _1 in our baseline model), which repre- 
sents residual turbulence due to breaking buoyancy waves (Lindzen 1981) and such. A list of 
prescribed/adjustable parameters is provided in Table 1. 

The remaining parameter in Equation 5 is the convective velocity scale from mixing-length 
theory: w* = K/L. Our baseline values for the turbulent mixing parameters just below the Jovian 
ammonia cloud are H = L = 20km, K = 2 x 10 8 cm 2 s _1 , and w* = 1.1ms -1 . 

To compute the vertical distributions of condensate and vapor, we proceed upwards from the 
sub-cloud conditions, requiring all excess vapor to condense and solving Equation 4 at each level. 
If we heuristically assume that q c /qt and L are constant in a cloud, the solution is an exponential 
decline of total mixing ratio with height above cloud base (where we define z = 0): 



V Qt V 

Note that by using the sub-cloud mixing ratio as a lower boundary condition, any moistening due 
to rain evaporating below the cloud base is ignored. 

Comparing our adaption of the Lunine et al. (1989) profile with our calculation using / ra j n = 1 
isolates the effect of reducing the mixing length due to atmospheric stability. The cumulative effect 
of the progressive reduction in mixing length due to the stability of the Voyager temperature profile 
above the cloud base is seen to result in a cold-trap in the lower stratosphere. Tripling / ra i n further 
reduces the cloud density and lowers the cold-trap to the tropopause; increasing it to 10 results in 
a cloud with less condensate than the Lewis (1969) model. 

We assume that / ra j n is independent of altitude. Yet specifying an appropriate value of / ra i n 
at the cloud base, let alone any vertical dependence, poses a significant challenge. For guidance, 
first we turn to in situ measurements and detailed simulations of terrestrial water clouds, and then 
consider constraints provided by values retrieved through remote sensing of Jovian ammonia clouds. 

For terrestrial stratocumulus clouds capping well-mixed planetary boundary layers, we find 
that / ra i n < 1 in the cloud deck and increases with distance below cloud top. An assortment of 



L = Fmax(A,r/T adiab ) 



(6) 




(7) 
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in situ measurements indicates that / ra j n increases with decreasing droplet concentrations (N), as 
fewer and hence larger droplets more efficiently produce drizzle and thereby decrease cloud water. 
For example, / ra i n ~ 0.2 for a case study over the North Sea, where N = 100 cm -3 (Nicholls, 1984), 
while in California stratocumulus it increased from 0.3 to 0.5 as N decreased from 40 cm~ 3 in clouds 
contaminated by ship exhaust, to 10 cm~ 3 in clean ambient air (from Figure 36 of Ackerman et 
al. 2000b). We note that in clean marine stratocumulus, reduced rain is observed (e.g., Taylor and 
Ackerman 2000) and predicted (e.g., Ackerman, Toon, and Hobbs 1993; Stevens et al. 1998) to 
result in deeper cloud layers, which is consistent with our simple model. However, any changes in 
cloud cover associated with precipitation changes are not represented in our model. 

We have also calculated profiles of / ra i n for deeper convection, using large-eddy simulations of 
trade cumulus clouds (Ackerman et al. 2000a), which are twice as deep as stratocumulus clouds 
(~1000 m compared to ~500 m). As in the stratocumulus clouds, we find that / ra i n < 1 in the 
stratiform anvils at cloud top. However, / ra j n is significantly enhanced throughout the bulk of 
the trade cumulus clouds, where it ranges from ~ 2 to 6. Vertical winds are more symmetrically 
distributed in stratocumulus, and hence the convective velocity scale is representative of the mean 
updraft velocity. The circulation is more skewed in trade cumulus, with narrow updrafts opposing 
the broad subsidence: for the Ackerman et al. (2000a) trade-cumulus simulations the mean vertical 
velocity in cloudy updrafts is three times the convective velocity scale. Such a skewness is consistent 
with an enhancement of the rain factor in deep convection. 

We recommend a more systematic analysis of / ra i n in terrestrial clouds. Perhaps more pertinent 
to the much deeper clouds expected in gas giants and brown dwarfs, we also recommend consid- 
eration of much deeper convection than considered here. For now we treat / ra i n as an adjustable 
parameter, leaning toward values > 1, as we expect the deep convection in substellar atmospheres 
(on the order of an atmospheric scale height and deeper) to more closely resemble cumuliform than 
stratiform convection in the terrestrial atmosphere. 

For observations of Jovian ammonia clouds, we first turn to the retrievals obtained from the 
Voyager IRIS (InfraRed Interferometer Spectrometer) instrument by Carlson, Lacis, & Rossow 
(1994), who considered latitudinal variations among zones and belts in the Jovian tropics. Al- 
though they did not provide profiles of condensed ammonia, we can compare our results to their 
ratios of condensate to atmospheric scale height (H p /H g in their notation). Carlson et al. (1994) 
retrieve ratios of 0.35 and 0.40 (± 0.10) for ammonia clouds in the Equatorial and Northern Trop- 
ical Zones, respectively, which are seen in Figure 3 to span a range of / ra i n values between ~ 1 
and 3. Retrievals of ammonia cloud properties for the Jovian tropics from ISO (Infrared Space 
Observatory) measurements by Brooke et al. (1998) indicate a scale height ratio of 0.3, which is 
consistent with our model results for / ra j n ~ 2. 

Were the ratio q c /qt fixed in our model, as assumed heuristically for Equation 7, the condensate 
height should vary as f~J n - However, as seen in Figure 3, the model dependence is not nearly that 
steep, and for / ra j n > 3 the dependence nearly vanishes. The dependence is moderated by a negative 
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feedback in which q c /qt decreases with increasing / ra i n - Our assumption of zero supersaturation 
within a cloud is equivalent to q c /qt = S — 1, where S is the potential supersaturation before 
condensation eliminates it. Increased / ra in offsets more of the turbulent mixing of vapor and 
condensate, thereby reducing the potential supersaturation and decreasing the ratio q c /qt- The 
negative feedback is thereby due to the reduction of q c /qt, which diminishes the dependence of the 
condensate scale height on / ra j n (Equation 7). 

The assumption that all vapor in excess of saturation condenses can be relaxed in our model 
by replacing Equation 1 with 

q c (z) = max(0, q v (z - Az) - {S ciouA + 1) q s (z)) (8) 

where S c i ou d is t ne supersaturation that persists after accounting for condensation. Allowing 
Scioud > represents conditions in which there is a significant barrier to the formation of cloud 
droplets, and/or condensation is too slow to effectively offset the supersaturation driven by cooling 
in updrafts. In shallow terrestrial water clouds, such as stratocumulus, neither of these condi- 
tions holds, as there are typically abundant condensation nuclei upon which droplets form at low 
supersaturations, and the concentration and diffusivity of water vapor are sufficient to allow con- 
densation to balance the modest dynamic forcing at low supersaturations (~10~ 3 ). However, in 
cirrus clouds that form directly from the vapor phase (as opposed to the freezing of water droplets 
from deep convection), there are typically few effective ice nuclei available, and hence the barrier to 
nucleation can result in supersaturation building to ~0.5 before ice crystals form, even in moderate 
updrafts. High supersaturations (~0.3) can be maintained after nucleation in cirrus clouds because 
the concentration and diffusivity of water vapor are greatly reduced at the cold temperatures of 
the upper troposphere (Jensen et al. 2000). 

As extreme cases, in Figure 2b the condensate profile for Sdoud = is compared to that for 
Scioud = 1- The barrier to condensation results in a lifting of the cloud base and enhanced lofting 
of vapor, but the altitude of the cold-trap is unchanged (recall that the temperature profile is fixed 
here). The greater lofting enhances condensate mass above ~0.4 bar, which tends to increase the 
column of condensate. However, the tendency is more than offset by the decrease in atmospheric 
mass density at the elevated cloud base, and hence the condensate column decreases from 52 to 
39 g m~ 2 in response to increasing S c i ou d from to 1. This 25% reduction of condensate column 
contrasts markedly with the effect of increasing / s from 0.01 to 1 in the Marley et al. (1999) model, 
in which the condensate column increases a hundredfold (Figure 2a). These opposite responses (to 
comparable changes in maximum relative humidities: from 100 to 200% for the present model, and 
100.01 to 200% for the Marley et al. model) arise from the distinct definitions of the supersaturation 
factor: Sdoud corresponds to the supersaturation after condensation is treated in the present model, 
whereas / s corresponds to the potential supersaturation before condensation is treated in the Marley 
et al. model. Marley et al.'s / s is more comparable to f^ n than to S c i ou d in the present model. 
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3.2. Droplet Size Distributions 



Size spectra of cloud particles in terrestrial condensation clouds are commonly observed as 
bimodal number distributions, with a condensational mode (of radius ~10 /xm in stratocumulus) 
resulting from condensational growth at modest supersaturations, and a precipitation mode at 
larger radii resulting from coalescence due to dispersion in sedimentation velocities. The bimodal 
structure is evident in droplet size distributions measured in stratocumulus clouds off the coast 
of California, as shown in Figure 4a. At the cloud base, the mean size of the condensation mode 
depends on factors such as the updraft velocity and the composition and size distribution of the 
condensation nuclei upon which droplets form. Above the cloud base, the mean size of the conden- 
sation mode increases with altitude in an updraft. 

The production of precipitation in terrestrial clouds is generally observed to increase with 
the mean size and spectral width of the condensation mode. The mean size of droplets in the 
precipitation mode 4 is found to increase with distance from cloud top in marine stratocumulus 
(e.g., Nicholls 1984) (the opposite tendency of the condensation mode), typically explained in 
microphysical terms as due to "fortunate" collector drops sweeping up smaller droplets and growing 
as they fall. Such a profile is consistent, as it must be, with the observed decrease in / ra i n with height 
above the cloud base. We note in passing that measured profiles of / ra i n are also consistent with our 
assumption that the size of precipitation particles decreases as convective velocity decreases above 
the cloud base (approaching the temperature inversion that caps the boundary layer), though we 
assume a uniform value of / ra j n in our model calculations. 

We make no attempt to model the complexity of cloud processes here, as such detailed com- 
putations are prohibitively demanding, and the parameter space of unknowns is overwhelming for 
the range of condensates expected in substellar atmospheres (including such basic issues as whether 
condensates are solid or liquid). Instead, we simply prescribe a single, broad lognormal size distri- 
bution of condensate particles at each level (Figure 4), thereby halving the number of parameters 
required for a bimodal distribution. The lognormal size distribution is given by 



where n is the number concentration of particles smaller than radius r, and the three parameters 
to be constrained appear on the right side: N is the total number concentration of particles, r g the 
geometric mean radius, and a g the geometric standard deviation. 

Marley et al. (1999) also prescribe a lognormal size distribution of condensate particles, in 
which <Tg is fixed at 1.5 and r g is determined from the particle sedimentation velocity corresponding 
to the convective velocity scale (w*). We choose to use a g as an adjustable parameter, and determine 



4 The maximum size of precipitation particles is set by their breakup due to hydrodynamic instability, which for 
terrestrial raindrops occurs at a radius of ~3 mm. 




(9) 
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N and r g through q c and / ra in- The rain factor is rigorously defined through 

_ Jo 00 vt (dm/dr) dr 

Jrain — l iu J 

where Vf is the particle sedimentation velocity (described in Appendix B), m is particle mass, and 
e is the ratio of condensate to atmospheric molecular weights. To close the system analytically we 
fit a power-law dependence for particle fallspeed about its value at Vf(r w ) = w* through 

v { = w*(r/r w ) a (11) 

where the exponent a is calculated from a fit to the fallspeeds between r w /a and r w when / ra i n > 1, 
and between r w and r w a otherwise (a is constrained to be > 1.1 for the fit). The power-law 
approximation allows Equation 10 to be expressed as 

_ Jo"" r 3+a (dn/dr) dr 

~ r« / °° r 3 (dn/dr) dr 1 j 

Integration of the lognormal distribution then leads to 



rw/ ra i n exp ( — In cj g ) (13) 



AT 3ep a g c / 9 2 \ 

Ar= 4^| eXP V "2 ln (14) 

where p p is the density of a condensed particle (see Appendix B). The parameter / ra i n can be 
interpreted in terms of microphysics by identifying the radius of mass- weighted sedimentation flux 
in Equation 12: 

_ f Jo°° r3+a (dn/dr) dr 
SCd V r 3 (dn / dr) dr 

which leads to / ra in = (r sc d/r w ) a . In Figure 46 the size distribution is weighted by precipitation, 
where it is seen that r se d < r w , consistent with / ra j n < 1. In contrast, for the large-eddy simulations 
of trade cumulus mentioned above, and as implied by the retrievals from the Jovian ammonia clouds, 
the droplets grows sufficiently large to satisfy r SC( j > r w . 

It is seen from Equations 12—14 that droplet sizes are decoupled from condensate mass (though 
both depend on / ra i n ); the condensate mass simply scales the distribution through N. Ignoring 
any vertical dependence of atmospheric stability (Equation 6), vertical variations in droplet sizes 
are due to the height dependence of convective velocity, leading to r w oc p]/ 3a , which yields a mild 
vertical dependence of approximately r w oc p]J A for our baseline Jovian ammonia cloud (for which 
a = 1.3, corresponding to a moderately turbulent sedimentation regime in which fallspeeds are 
reduced from those in viscous flow). 
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Of greater interest than the mild vertical dependence are the sensitivities of droplet sizes on 
/ ra i n and <T g at the base of our Jovian ammonia cloud, where r w = 35 fim. For a given value of <j g , 
more efficient rain implies larger droplets: r g oc f]^- For our baseline droplet distribution (with 
Ug = 2), a value of / ra i n = 3 leads to r g = 14 fim (compare to ~ 10 fim for shallow terrestrial 
clouds), while / ra j n = 5 yields r g = 20 ^m. 

For efficient precipitation (/ ra i n > 1)> reducing the width of the size distribution requires 
an increase in r g , as the narrower distribution is centered on a larger radius. For example, a 
monodisperse distribution (cr g = 1) with / ra i n = 3 results in r g — 77 /xm. Such narrow distributions 
of large particles could result from condensational growth at high supersaturations, reminiscent of 
methane "rain without clouds" suggested by Toon et al. (1988) for Titan's atmosphere. 

In the complete model, we calculate spectrally-resolved profiles of condensate opacity by in- 
tegrating the scattering and absorption coefficients (from Mie calculations) over the particle size 
distributions. Here we simply present opacities for geometric scatterers, which for a model layer of 
thickness Az is given by 

Ar= 3e Ml ^ 
2 p P r cS 

where the effective (area-weighted) droplet radius is evaluated from the lognormal size distribution: 

r cS = TV/rYin exp (~^y^ In 2 ^ (17) 

The effective radius can be greater or less than r w , depending on the combination of / ra i n , a, and 
<r g . At the base of our baseline Jovian ammonia cloud r c g = 46 fim, which is 11 /im greater than 
r w . 

The computation of cloud optical depth depends on vertical grid resolution due to the exponen- 
tial temperature dependence of saturation vapor pressures. To reduce such resolution dependence 
we progressively subdivide each model layer until its optical depth converges to 1% precision. For 
the Jovian ammonia cloud our calculations converge at a minimum sub-layer thickness of ~ 30 m. 

Equations 16 and 17 show that increased precipitation reduces opacity not only by decreasing 
q c but also by increasing r e ff . For heuristic purposes the column optical depth can be estimated from 
the cloud base properties by ignoring any height dependence of the mixing length and atmospheric 
scale height and assuming q c = q t within the cloud, in which case 

_ 3 epgbciow /,o\ 
2 5 r eff (l + / rai „) 

where p and r e g are atmospheric pressure and droplet effective radius at the cloud base. Note that 
r is more than linearly dependent on % e i OW ) due to its dependence on the cloud base pressure. For 
the following comparisons with observations of the Jovian ammonia cloud, we sum optical depths 
from Equation 16 over the model layers and compaute a cloud average r e g- from the cloud optical 
depth and vertical column of condensate. 
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4. Comparisons with Observations 

West et al. (1986) attempt to reconcile among a vast array of observations of the Jovian 
ammonia cloud, and conclude that its optical depth at visible wavelengths is between ~ 2 and 10, 
comprised of a population of small particles (r ~ 1 jum) reaching the tropopause in both belts and 
zones, underlain in zones by a population of larger particles (r ~ 3 — 100 fj,m) concentrated near 
the cloud base. 

More recently, Banfield et al. (1998), who make no attempt to retrieve particle sizes, assume 
a particle effective radius of 0.2 /xm and retrieve optical depths from Galileo imaging data at a 
wavelength of 0.756 that cluster in the range ~ 1 — 4.5 (their Figure 9). For comparison with 
the optical depths of ours and West et al.'s (at a mid-visible wavelength of 0.55 /Ltm), we scale 
Banfield et al.'s r by the ratio of extinction efficiencies at 0.55 to 0.756 /jm for 0.2 /im particles 
(using their refractive index of 1.4), resulting in a mid-visible r range of ~ 2 — 10. Hence, the West 
et al. (1986) and Banfield et al. (1998) results are effectively identical, and are hereafter lumped 
together as "West et al." 

Recall from Figure 3 that the condensate scale height retrievals of Carlson et al. (1994) are 
consistent with our baseline model for / ra i n between ~ 1 and 3. This entire / ra i n range for our 
baseline model (in which cr g = 2) is seen in Figure 5a to overlap with the data reported by West 
et al. (1986). The overlap corresponds to r c g ranging from ~ 18 to 45 ^m. A narrower size 
distribution results in larger droplet effective radii and therefore smaller optical depths; for the 
monodisperse case the data mutually overlap at r e ff ~ 34 to 60 /im. 

Carlson et al. (1994) use Voyager IRIS spectra to retrieve cloud optical depths 5 and droplet 
sizes for the Northern Equatorial Belt 'hot-spots' (note that the infrared data are insensitive to 
the submicron particles reported by West et al. 1986). Their best-fit particle distribution is a 
mixture of small (r e ff = 3 ^m, r = 0.16) and larger (r e fr = 100 /im, r = 0.38) particles, resulting 
in a combined r c g = 72 /im and r = 0.54. However, hot spots are anomalous features of reduced 
cloudiness associated with pronounced dynamical forcings (e.g., Showman & Ingersoll 1998), and 
our one-dimensional model is intended to represent horizontally averaged conditions, which would 
seem more comparable to the Equatorial and Northern Tropical Zones, where Carlson et al. retrieve 
r of ~ 1.2 and 2, respectively. Carlson et al. do not present separate retrievals of droplet sizes 
for the zones, but do state that their optical depths are dominated by larger particles, which we 
interpret as an observed range of 70 fj,m < r c g < 100 ^m. 

Given that our model results overlap with the entire range of the West et al. (1986) optical 
depths, which in turn do not overlap with those of Carlson et al. (1994), it should not be surprising 
that there is no mutual overlap between our baseline model results and the retrievals of r and r c g 



5 Carlson et al. (1994) report cloud optical depths for an assumed extinction efficiency of 1, whereas by treating 
cloud particles as geometric scatterers, we assume an extinction efficiency of 2. Hence, we multiply their optical 
depths by 2 for comparison with ours. 
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by Carlson et al. (1994). However, our model results with / ra i„ = 3 just overlap with the Carlson et 
al. (1994) retrievals if a monodisperse size distribution is assumed, which suggests that the Carlson 
et al. data are dominated by particles in the precipitation mode. Alternatively, multiplying our 
eddy diffusion coefficients by a factor of 3 also results in larger particles and therefore reduced 
optical depths, which leads to overlap with the Carlson et al. (1994) retrievals for / rain ~ 2 (Figure 
56). 

The above comparisons show that reasonable choices of model parameters produce agreement 
with observations of tropical Jovian ammonia clouds. The few unknown model parameters in our 
simple model are not uniquely constrained, befitting the incompleteness and uncertainty in the ob- 
servations and the ambiguity in comparing a model for globally averaged clouds with measurements 
that resolve large-scale horizontal variability in the clouds. To compare our model with observa- 
tions averaged over a wider area we next consider measurements of the Jovian tropics obtained 
from Earth orbit. 

Brooke et al. (1998) use a 3-fim ISO spectrum to retrieve microphysical properties of the 
Jovian ammonia cloud and find best fits for two possibilities: first, a monomodal distribution of 
10-/mi ammonia particles with a visible optical depth of 1.1; and second, a bimodal distribution of 
1 and 10-//m ammonia particles with an optical depth of 1.3 (equivalent to an effective radius of 7 
(im). Both fits include an additional optical depth of 0.1 from grey particles. In their analysis, a fit 
for 10-/xm particles indicates a superior fit in comparison to those for 1 and 30 fim particles, which 
we intepret as allowing a size range of 5 to 20 /im. Recalling from the discusson of Figure 3 that the 
Brooke et al. (1998) condensate scale height (a single value with no uncertainty) is consistent with 
/rain = 2. Figure 56 indicates that mutual consistency with the effective radii retrieved by Brooke et 
al. (1998) requires increasing the width of our size distrubutions and/or reducing our eddy diffusion 
coefficients. Simultaneously matching the Brooke et al. (1998) optical depths requires a substantial 
reduction of the sub-cloud ammonia abundance. However, the baseline ammonia abundance we 
use is reported by Brooke et al. (1998) (as a single value with no uncertainty) to best fit their 
3-//m spectrum. Hence, our calculations are evidently inconsistent with that baseline ammonia 
abundance and the combination of small particles and small optical depths retrieved by Brooke et 
al. (1998). 

We have already noted the shortcoming that our model excludes the possibility of horizontal 
variability. For the case of modeling emitted radiative flux, this simplification will of course result 
in an underestimate at some wavelengths, since any flux leaking out through the clearings between 
patchy clouds is not treated. Horizontal variability is also ignored in the retrievals of Brooke et al. 
(1998), which results in their underestimating cloud optical depth due to a plane-parallel albedo bias 
(e.g., Calahan et al. 1994): the area- weighted albedo of a cloud deck calculated from a single column 
with an optical depth r is always greater than the albedo averaged over a variety of columns with 
the same area-weighted average r. Inverted for the Brooke et al. (1994) retrievals of optical depth 
from reflected spectra, this bias indicates that the optical depth from an area-weighted average 
radiance underestimates the area-weighted optical depth. Hence, the actual area-weighted optical 
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depth is greater than reported by Brooke et al. (1994) and closer to the optical depths calculated 
by our model. We are unable to quantify the magnitude of this error without reproducing their 
retrievals, which is beyond the scope of this study. 

5. Applications to Substellar Atmospheres 

5.1. Vertical Cloud Structure 

We use our baseline model (Table 1) to calculate profiles of condensed water, silicate (as 
enstatite, MgSiOa), and iron in theoretical atmospheres of brown dwarfs and a giant planet (see 
Figure 6 for gravities and effective temperatures). The temperature profiles are calculated for 
cloud-free conditions (Marley 2000). The L-dwarf-like atmosphere (Figure 6a) is too warm for 
water to condense. Between 1 and 10 bars, within the convective region, the silicate and iron 
clouds are seen to overlap, suggesting the possibility of microphysical interactions between them, 
which are ignored by our model. The silicate particles in this case are about twice as large as the 
iron particles because the assumed density of a silicate particle is about half that of an iron particle 
(Appendix B). 

Although the temperature profiles in Figure 6 were not calculated self-consistently to include 
the effects of clouds, it is clear that for objects near T c g ~ 1500 K, clouds are an important 
opacity source. The silicate and iron clouds in the L-dwarf-like model (Figure 6a) appear in the 
visible atmosphere and therefore play an important role in controlling opacity and the temperature 
structure of the atmosphere. Nevertheless these clouds are confined to a relatively thin cloud 
deck, which does not reach the upper regions of the atmosphere as do condensates in the well- 
mixed profiles also shown. The cloud particles are also fairly large (r e ff ~ 40 — 80 fj,m) and will 
have a substantially different spectral opacity than smaller particles. The Lyon group (Chabrier 
et al. 2000) employs an "astrophysical dust" size distribution of sub-micron particles to model 
dust opacity in such atmospheres. The cloud model presented here, with larger particles confined 
to a discrete cloud deck, represents a substantial departure from the previous work. In a future 
publication we will discuss the spectral and color properties of atmospheres with these new cloud 
models. 

Figure 66 presents the cloud model applied to a T-dwarf-like atmosphere with T e $ = 900 K, 
which is again too warm for water to condense. (No iron cloud is shown in Figures 66 and c 
because the cloud base is below the bottom of the model domain.) Although the silicate cloud and 
the omitted iron cloud may be important to the atmospheric temperature structure, they no longer 
represent significant opacity sources to an observer. 

The changing role of cloud opacity with effective temperature is more clearly shown in Figure 
7, which illustrates the brightness temperature spectra of several radiative-equilibrium models for 
brown dwarf atmospheres as well as the atmospheric temperature range over which most of the 
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cloud opacity is found. In a model with T e g = 1800 K (Figure 7 a) the silicate cloud deck forms in 
the model stratosphere and is relatively thin. Comparison of the solid and dotted curves, which 
respectively include and exclude silicate and iron cloud opacities, shows little difference between 
the two cases. Since the cloud optical depth is only a few tenths, flux is efficiently transported 
from levels deeper than the base of the cloud. The iron cloud (not shown) adds a few more tenths 
of optical depth. Hence, the clear and cloudy models are very similar. Such a model would be 
appropriate for an early type L-dwarf. 

Figure 7b shows the results for a cooler atmosphere, with T e g = 1400 K, appropriate for a late 
L-dwarf (Kirkpatrick et al. 1999; Stephens et al. 2001). Here the cloud is much more optically 
thick, and flux originates no deeper than the middle of the cloud layer. In the clear atmosphere, 
flux originates from deeper, hotter levels. As a result the band depths of the cloudy model are 
shallower, a result shown for dusty M-dwarfs by Jones and Tsuji (1997). Since flux is conserved for 
T e g fixed, the cloudy model emits more flux beyond about 2 ^m than the clear model. Regions of 
strong molecular absorption, including the depths of the water bands shortwards of 2 /xm as well as 
most of the 2 to 5 fim region evidence higher brightness temperatures in the cloudy case since the 
atmosphere above the cloud must warm to produce the same total emitted flux as in the cloud free 
calculation. Regions of stronger molecular opacities are sensitive to these warmer temperatures 
higher in the atmosphere. 

In Figure 7c, for which T e g- = 900 K the silicate cloud forms well below the region in which 
most flux originates and again the clear and cloudy models are similar. However in the regions in 
which the molecular opacity is lowest, near 1.1 and 1.3 fxm, flux originates from deeper regions in 
the clear atmosphere than for the cloudy case. As a result the peak-to-trough variation in emitted 
flux is again somewhat smaller for the cloudy model. We note that clear atmosphere models for 
T-dwarfs like G1229 B and GD165B typically over predict the water band depths (e.g., Marley et 
al. 1996; Allard et al. 1996; Tsuji et al. 1996; Saumon et al. 2000; Geballe et al. 2001), and 
suggest that the attenuation of flux by the top of the silicate cloud deck may be responsible for 
this effect. 

Notably, these model atmospheres illustrate the origin of the curious change in infrared colors 
of the L- and T-dwarfs (e.g., Kirkpatrick et al. 1999; Martin et al. 1999; Fan et al. 2000). The 
cloud-free cases shown in the figure monotonically vary in J — K from 1.39 to -0.17 from warmest to 
coolest. In contrast the cloudy models initially become redder with falling T c g (moving from 1.6 to 
1.7) before they move to the blue ( J — K = 0.38 for the case in Figure c) when the silicate and iron 
clouds begin to disappear below optical depth unity in the gas. Thus our precipitating condensation 
cloud model qualitatively reproduces the color variation of the L- and T-dwarfs, consistent with 
previous arguements based on interpretation of spectra (Allard et al. 1996; Marley et al. 1996; 
Tsuji et al. 1996). In contrast, the pure chemical equilibrium model (Chabrier et al. 2000) predicts 
the presence of substantial dust opacity well to the top of the atmosphere, which is clearly excluded 
by the data. A more complete treatment of color changes will be given in a future study. 



-19- 



For the cooler atmosphere representative of a cool extrasolar giant planet (Figure 6c), water 
condenses in the radiative region, in essence a stratospheric cloud. The ice particles are seen to 
be larger than the silicate particles, chiefly due to the lesser densities assumed for the individual 
ice particles. The reduced gravity in the atmosphere of the less massive extrasolar giant planet 
requires larger silicate particles to match the mean sedimentation velocity than does the more 
massive T-dwarf. Note again that the well-mixed assumption produces a profoundly different 
vertical structure, in which the silicate cloud is so deep that it significantly overlaps the water 
cloud. 

As did Mar ley et al. (1999), Sudarsky et al. (2000) have computed water cloud profiles in order 
to estimate extrasolar giant planet albedos. Sudarsky et al. (2000) essentially assume / ra i n = 
and limit the cloud to be no more than 1 scale height thick. Such a model would be similar to the 
well-mixed water cloud in Figure 6c with a flat cloud top at ~ 4 x 10~ 3 bars. 

The emergence of water clouds in substellar atmospheres with T e g below about 500 K will 
reshape the vertical temperature profile and emergent spectra of these objects. Preliminary models 
computed with this cloud profile suggest that such cool objects will again move to the red in J — K 
after the blueward excursion caused by the sinking of the silicate cloud below the visible atmosphere 
and the emergence of CH4 as a dominant opacity source in K band. Hence the near-IR colors of 
very cool objects computed from cloud- free atmosphere models (e.g., Burrows et al. 1997) are likely 
to differ substantially from actual objects. 

5.2. Non-uniform Clouds 

Variable brightness in / band has been detected for some L-dwarfs by Bailer-Jones & Mundt 
(2001), who attribute the variability to evolution of dust clouds. They find some evidence that 
variability may be more common in later-type L-dwarfs. Although we do not model horizontally 
variable clouds, these observations are consistent with the model presented here. As clouds form in 
progressively cooler objects they become more optically thick and form deeper within the convective 
region of the atmosphere. Thus global scale tropospheric weather patterns, as seen on Jupiter 
and predicted for brown dwarfs (Schubert & Zhang 2000), can more easily produce photometric 
variability since the turbulent motions are greater, making local clearings more likely, and enhancing 
the potential contrast between clear and cloudy air. Indeed the great red spot of Jupiter produces 
a photometric signal in both reflected sunlight and emitted thermal radiation (Gelino & Marley 
2000). 

Horizontally varying silicate clouds, even if not of the appropriate scale to produce a varying 
photometric signal, may play an important role in the transition from the dusty L-dwarfs to the 
relatively cloud-free T-dwarfs. The change in J — K color from the latest red L-dwarfs ( J — K ~ 2) 
to the blue T-dwarfs (J — K ~ 0) is quite abrupt. Four L8 dwarfs with known or estimated 
absolute magnitudes are only 1 magnitude brighter in J band (Reid et al. 2001) than G1229B. 
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Reid et al. argue that this implies the L8 dwarfs are only about 250 K warmer than G1229B. 
Even with the silicate cloud deck forming at progressively deeper levels with falling T e g , it may be 
difficult to account for such rapid color variation. In fact the rapid transition may be a signature 
of horizontally varying clouds. Once tropspheric convective patterns begin to produce substantial 
horizontal variability, the flux from the more cloud-free regions will begin to dominate the total 
emitted flux, even if large fractions of the object are still cloudy. For example, Jupiter's 5-fim flux 
is dominated by the relatively cloud-free 'hot-spots' (Westphal et al. 1974) that typically cover 
about 1% of the surface area of the planet (Orton et al. 1996). Thus the apparent rapid change 
from cloudy L-dwarfs to clear T-dwarfs may be due to a gradual change in cloud coverage in the 
visible atmosphere, with the larger flux from the clear regions quickly dominating. 



We have developed a simple cloud model for substellar atmospheres that includes precipitation 
by condensate particles larger than that set by the convective velocity scale, which permits us 
to reproduce the properties retrieved from Jovian ammonia clouds. Effective precipitation also 
produces cloud profiles in theoretical brown dwarf and extrasolar giant planet atmospheres that 
are broadly consistent with observations. 

As in the solar system, real clouds in the atmospheres of substellar objects will likely be neither 
uniform nor homogeneous, however we hope that this model will provide a framework for evaluating 
the globally-averaged role such clouds play in controlling the thermal radiative transfer and spectra 
of brown dwarfs and extrasolar giant planets. 
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We also thank Robert West and Kevin Zahnle for providing helpful comments on the manuscript. 
M.S.M. acknowledges support from NASA grants NAG58919 and NAG59273 as well as NSF grants 
AST 9624878 and AST 0086288. 



For the saturation vapor pressure of ammonia (e s , in dyne cm 2 ) we fit the measurements 
tabulated in the CRC handbook (Weast, 1971) with 



where the temperature is in K. 

For the vapor pressure of water we use the expressions of Buck (1981), over ice for T < 273.16 K, 



6. Summary 



A. Saturation Vapor Pressures 
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and over liquid water at greater temperatures: 

e s (H 2 0,ice) = 6111.5exp (A2a) 

e s (H 2 0, liquid) = 6112.1 exp( 18 ' 72 ^;^ 8 / 7 227 - 3 ) (A2b) 

where T c is the temperature in degrees Celsius. These expressions are unsuitable at T > 1048 K, 
leading to vapor pressures that decrease with with increasing temperatures. Hence, at greater 
temperatures we simply fix e s (H20) = 6 x 10 8 dyne cm" 2 , which is its value at T = 1048 K. 

The vapor pressures for iron and enstatite are taken from Barshay and Lewis (1976). For iron 
below and above its melting point of 1800 K, we use, respectively 

/ 47664 \ 

e s (Fe, solid) = exp M5.71 — J (A3a) 

( 37120 \ 
9.86 — J (A3b) 



and for enstatite we use 

( _ 5. 

— ) 



e s (MgSi0 3 ) = exp ( 25.37 - (A4) 



B. Sedimentation Velocities 

Droplet terminal fallspeeds are calculated by first assuming viscous flow around spheres cor- 
rected for gas kinetic effects: 

2 Pgr 2 Ap 

vt = Bl 

9 7] 

where g is the gravitational acceleration, r is the droplet radius, and Ap = p p — p& is the differ- 
ence between the densities of the condensate and the atmosphere. The Cunningham slip factor, 
(3 = (1 + 1.267V Kn ), accounts for gas kinetic effects, in which the Knudsen number (-/Vkii) is 
the ratio of the molecular mean free path to the droplet radius. The dynamic viscosity of the 
atmosphere is given by Rosner (2000): 

5 y^fcjr (k B T/ef w 
71 =16^1* L22 - (B2) 

where d is the molecular diameter and e is the depth of the Lennard-Jones potential well for the 
atmosphere (2.827 x 10~ 8 cm and 59.7&b K, respectively, for H 2 ) and fee is the Boltzmann constant. 

For turbulent flow, at Reynolds numbers (N^ e = 2rp a Vf/r]) between 1 and 1000, we use a 
standard trick to solve the drag problem. Noting that CdN^ c = 32p a gr 3 Ap/3i] 2 is independent of 
fall velocity, we fit y = log(iVRe) as a function of x = log(CdiV^ c ) to the following data: at N-& e 
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= 1 we assume viscous flow, with = 24; for intermediate Reynolds numbers we use the data 
for rigid spheres from Table 10-1 of Pruppacher k Klett (1978); and at 7Vr 6 = 1000 we assume 
an asymptote of Cd = 0.45. This asymptote is appropriate to moderately oblate spheroids (Figure 
10-36 in Pruppacher k Klett 1978), which are more appropriate to unknown condensates than the 
extreme case of smooth spheres. Our fit to the data is y = 0.8x — O.Olx 2 , which allows us to evaluate 
the droplet terminal fall velocity from Nn c . 

At Reynolds numbers > 1000 we assume the drag coefficient is fixed at its asymptotic value 
(Cd = 0.45), which leads to 



We assume rigid particles and thereby ignore breakup, for instance by liquid droplets due to 
hydrodynamic instability. For the density of ammonia ice particles, we use 0.84 g cm" 3 (Manzhelii 
k Tolkachev 1964); for water we use 0.93 g cm" 3 (corresponding to ice at a temperature of 200 K, 
using Equation 4-17 from Pruppacher k Klett 1978); and for enstatite and iron we use 3.2 and 7.9 
g cm" 3 , respectively (Table 1.18 of Lodders k Fegley, 1998). 
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Table 1. Adjustable model parameters. 



Parameter Baseline value Description 



/ ra i n 3 Ratio of mass-weighted sedimentation velocity to convective velocity scale 

^"min 10 5 cm 2 s~ 1 Minimum value of eddy diffusion coefficient 

A O.f Minimum ratio of turbulent mixing length to atmospheric scale height 

Scloud Supcrsaturation that persists after accounting for condensation 

Og 2 Geometric standard deviation in lognormal size distributions of condensates 
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Fig. 1. — Vertical profiles of mole fraction (mixing ratio by volume) of condensed ammonia (q c ) 
from present model of Jovian ammonia cloud with different values / ra in 5 and from our adaptations 
of other models as labeled. The vertical coordinate is atmospheric pressure. The dotted line is the 
temperature profile. The kinks in the condensate profiles are due to ripples in the temperature 
profile. 
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Fig. 2. — Vertical profiles of condensed ammonia (as in Figure 1) from (a) the model of Marley 
et al. (1999) for two values of / s (the potential supersaturation prior to condensation), and from 
(b) our baseline model for two comparable values of 5 c i ou d (the supersaturation persisting after 
condensation) . 
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Fig. 3. — The ratio of the condensed ammonia scale height to the atmospheric pressure scale height 
(H c and H, respectively) as a function of / ra in- The model condensate scale height is calculated 
from the altitude of the peak opacity (for geometric scatterers) and the altitude at which the opacity 
falls to exp(-l) of its peak. The grey region depicts the range of scale height ratios retrieved from 
the Jovian Equatorial and Northern Tropical Zones by Carlson et al. (1994). 
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Fig. 4. — Measured size distributions (squares) of (a) cloud droplet number concentration (n, in 
cm~ 3 ) and (b) precipitation flux (P, in mm d _1 ) measured in California stratocumulus and averaged 
over 12 km of flight at 200 m altitude, ~100 m below cloud top (from Figure 3b of Ackerman et al. 
2000b). Solid curve is a bimodal fit to the measurements, the sum of two lognormal distributions 
(dotted curves). Dashed curve is a monomodal lognormal fit to the measurements, in which N = 
40 cm -3 , r g = 7 fim, and a g = 1.8 (symbols defined in text). The values of r g , r e s, and r se( j are 
tied to the monomodal fit to the measurements, while r w is calculated from the convective velocity, 
measured to be 0.33 m s _1 (I. Brooks 1995, private communication). 
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r eff (m) 

Fig. 5. — Condensate optical depth (for geometric scatterers) plotted as a function of particle 
effective radius (r e s) for the Jovian ammonia cloud. The darkest grey region corresponds to the 
range of observations given by West et al. (1986) (equivalent to Banfield et al. 1998, as described 
in the text), the medium grey region corresponds to the retrievals from Voyager IRIS 5 — 45 ^m 
observations by Carlson et al. (1994), and the light grey region corresponds to retrievals from ISO 
3-/UH1 observations by Brooke et al. (1998). The lines correspond to variation of model results 
as /rain ranges from 0.1 to 10, with results at / ra j n of 1, 3, and 5 marked by symbols as denoted 
in (a). The particle effective radii correspond to an opacity- weighted averages. The solid lines 
correspond to baseline values of (a) the width of the log-normal particle size distribution (a g ) 
and (b) the eddy diffusion coefficient (K). The dotted line in (a) corresponds to a monodisperse 
particle size distribution; the dotted and dashed lines in (b) correspond respectively to dividing and 
multiplying K by a factor of 3. Filled- in squares denote model results using the complete baseline 
set of parameters. 
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Fig. 6. — Profiles of temperature (dotted curves) and condensate mass concentration from baseline 
model (p c , solid curves) in theoretical atmospheres of a (a) L-dwarf, (b) T-dwarf, and (c) extrasolar 
giant planet. Droplet effective radii at cloud base are shown as horizontal bars. Well-mixed clouds 
are shown as dashed curves. The theoretical temperature profiles are calculated for cloud-free 
conditions. 
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Fig. 7. — Brightness temperature as a function of wavelength for atmosphere models calculated 
self-consistently (Marley et al. 2001) to include (solid) or exclude (dotted) silicate and iron clouds. 
Brightness temperature increases downward to indicate increasing depth in the atmosphere. Clouds 
are calculated using the baseline model parameters. Solid straight line indicates base of silicate 
cloud, dashed line denotes level in atmosphere at which column extinction optical depth reaches 
0.1, and shading depicts the decrease in cloud extinction with altitude. Since cloud particle radius 
exceeds 10 /xm in these models, the Mie extinction efficiency is not a strong function of wavelength 
over the range shown. Shown are models characteristic of (a) an early- type L-dwarf with T e g- = 
1800 K, (b) a late L with T cff = 1400 K, and (c) a T-dwarf with T cS = 900 K. All atmosphere 
models are for solar composition and gravitational acceleration of 1000 m s -2 , roughly appropriate 
for a 30 Jupiter-mass brown dwarf. Sodium and potassium lines, calculated using the theory of 
Burrows, Marley, and Sharp (2000), are prominent in the optical. 



